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Wind  energy  has  many  advantages,  it  does  not  pollute  and  it  is  an 
inexhaustible  source.  However,  the  cost  of  this  energy  is  still  too  high  to 
compete  with  traditional  fossil  sources.  The  yield  of  a  wind  turbine  depends 
on  three  parameters:  the  power  of  the  wind,  the  turbine  power  curve  and  the 
ability  of  the  generator  to  respond  to  fluctuations  in  the  wind.  This  article 
presented  the  MPPT  of  a  wind  turbine  system  equipped  with  an 
asynchronous  generator  has  dual  power  under  Matlab  Simulink  program,  in 
the  first  time  we  simulated  all  the  conversion  chain  with  complete  model  of 
DFIG  and  vector  control  in  second  step  then  applied  the  extracted  maximum 
power  MPPT  strategists,  this  command  is  effective  and  has  several 
advantages  it  offered  to  kept  the  maximum  power  delivered  to  network 
despite  all  the  parameter  is  change. 
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1.  INTRODUCTION 

Electrical  energy  is  crucial  to  any  socio-economic  development  [1].  The  demand  for  electricity  is 
very  important  to  this  dilemma,  it  is  necessary  to  appeal  to  new  energy  sources  that  will  be  without 
consequence  for  man  and  environment.  Wind  energy  represents  a  sizeable  potential  for  bearing  damping 
demand  increasingly  rampant,  after  centuries  of  evolution  and  further  research  in  recent  decades  and  some 
wind  power  projects  developed  by  major  central  wind  turbines  provide  electricity  in  parts  of  the  world  at  a 
competitive  price  than  the  energy  produced  by  conventional  plants  [2] . 

Today,  the  development  and  proliferation  of  wind  turbines  have  led  researchers  in  Electrical 
Engineering  to  conduct  investigations  in  order  to  improve  the  efficiency  of  electromechanical  conversion  and 
quality  of  the  energy  supplied  [3].  It  is  in  this  context  that  this  paper  describes  that  we  present  a  study  on  the 
use  of  asynchronous  machines  has  dual  power  in  a  wind  system.  Initially  modeling  and  conversion  system 
simulation  (turbine  and  asynchronous  generator  dual  feed) 

Then  the  vector  control  in  stator  active  and  reactive  powers  are  proposed  in  a  last  step  And 
improved  the  system  study  and  extract  maximum  power  presents  MPPT  strategists  and  the  simulation  results. 
The  objective  of  this  paper  is  to  develop  a  novel  combined  MPPT -pitch  angle  control  system  of  a  variable 
speed  wind  turbine. 
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2.  MODELLING  WIND  POWER  SYSTEM 

2.1.  Model  of  the  wind  turbine 

The  wind  turbine  collects  the  kinetic  energy  of  the  wind  and  converts  it  into  a  torque,  which  turns 
the  blades  of  the  rotor  [3].  In  our  study,  this  system  used  a  wind  turbine  based  on  DFIG,  which  supplies  to 
generate  the  torque  required  to  the  load  in  Figure  1 . 


Gearbox 


Generator 


Figure  1 .  Schematic  of  the  wind  turbine 


2.1.1.  Wind  modeling 

In  our  case,  the  wind  speed  will  be  modeled  as  a  sum  of  several  harmonics  [5]: 

V(t )  =  10  +  0.2  sin(0.1047t)  +  2  sin(0.2665t)  +  sin(1.2930t)  -  0.2  sin(3.6645t)  (1) 


2.1.2.  Aaerodynamic  power 

Pt  =  i.Cpap).p.S.v3  (2) 

Which  0  is  defined  by:  A  =  (3) 

The  coefficient  of  power  Cp  represents  the  aerodynamic  efficiency  of  the  wind  turbine.  We  will  use 
an  approximate  expression  of  the  power  coefficient  given  by  [6] : 

Cp(A,/?)  =  (0.35  -  0.00167)).  (/?  -  2).  sin  “  0.000184.  (A  -  3).  (/?  -  2)  (4) 

The  simulation  of  the  power  coefficient  is  shown  in  Figure  2: 


Figure  2.  Cp  (A.,13)  Characteristics  for  various  values  of  pitch  angle  fi 


The  aerodynamic  torque  is  expressed  as  follows: 


ct  =  o  =  ttP-  r3- y2- CP(A,  (3) 


(5) 
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For  the  mechanical  part  can  thus  lead  to  a  more  simple  mechanical  model  [7  and  14]  (Figure  3).  The 
dynamic  behavior  of  the  generator  can  be  represented  by  the  following  equation: 


-T  -  f  w 

em  Jm 


(6) 


Where  Jt  is  the  rotational  moment  of  inertia  of  the  rotor  and  the  generator  kg.m2,  co  is  the  angular 
velocity  of  the  rotor  in  red  /  s,  Te  is  the  mechanical  torque  applied  to  the  alternator  shaft  in  Nm,  Tem  is  the 
electromagnetic  torque  developed  by  the  generator  in  Nm  and/m  is  the  viscous  friction  coefficient  in  Nm. 


Figure  3.  Mechanical  model  of  turbine 


2.2.  DFIG  Modeling 

We  modeled  the  DFIG  with  the  implementation  of  the  transformation  of  the  park  following 
Repository  related  rotating  field  This  repository  is  called  system  of  axes  (X,  Y),  it  rotates  with  the  speed  of 
the  electromagnetic  field  is  obtained  [9] : 


The  stator: 


The  rotor: 


(  —  D  ■  .  d(Pds 

uds  —  ttslds  H  ~  ^s^Pqs 

_  ■  d(pQS 

<  uqs  —  P^qs  “I  ^  1“  ^s^Pds 


Uos  P-s 


dffios 

“dT 


udr  —  Pr^dr  T  ^  (Ci)s  U)r)(Pqr 
uqr  =  Priqr  ^  (ws  —  wr)(Pdr 


dtpor 

dt 


(7) 


(8) 


The  equations  of  flow  and  after  the  simplification  are: 


fcpds  —  fs^ds  T  Midr 

j  ^Pqs  =  hs^qs  T  Miqr  (9) 

V  ^Pos  —  hs^os 


Also  for  the  rotor  and  in  the  same  way  are: 


<Pdr  —  Lri(ir  T  Mids 

( Pqr  ~  Lriqr  T  MiqS  (10) 

(Por  —  Lrior 
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Doubly-Fed  Induction  Generator  model  of  the  landmark  Park  in  the  form  of  state: 

~dt~  ~  LS(J  ^ds  ~  3  J^—idr  F  (oJsCT  —  PQr((J  —  1  ))hsiqS  —  J^^dr  F  Pf2rMigr) 

fit  ~  ~L~G  ^  X  (^s<7  —  Pf^rC^  —  I ))l"sGs  —  ~j^~^dr  F  Pf2rMiqr) 

didr  1  .  MPr  .  <r-l  .  M  a)sa  -  P£lr  . 

=  T~a^Udr  ~  ^r^dr  ^  £  ^ds  ^  iv/  PPLrLsLriqs  —  J^~uqs  3  ^  iqr) 

~dt~  ~  Lr( 7  (U<lr  ~  3  F  (&)s0'  +  Pf2r)LrGr  —  ~j^UqS  +  PflrMGs) 

_  _  3PM  ,  .  _  .  n  Cr  fr 

V  dt  ~  2Lr  Wdslw  Vcisldr)  J-  J-  Ur 

M2 

Where:  a  =  1 - is  the  dispersion  coefficient. 

LsLr 

The  expression  of  the  electromagnetic  torque  of  the  DFIG  in  reference  Park  [10].  The  general  form 
of  electromagnetic  torque  is: 

Qm  =  -^^{SPdriqs  ~  ^Pqr^ds)  =  (SPds^qr  ~  ^qs^dr)  (H) 


2.3.  Power  Modeling  DFIG 

2.3.1.  Modeling  of  phase  rectifier  diodes 

If  we  neglect  the  effect  of  encroachment,  the  rectifier  output  voltage  will  be  defined  as  following: 
Vred(t)  =  Max [Va (t) .  Vb (t) .  Vc (t) ]  -  Min [Va (t) .  Vb (t) .  Vc (t) ]  (12) 

2.3.2.  Modeling  of  DC  bus 

To  reduce  the  ripple  of  the  voltage  source  adding  a  low  pass  filter  LC  Their  operation  is  governed 
by  the  following  equations: 


(dId  _  Vred-VDC 
dt  Lf 

dVDC  _  Id-h 

dt  Cf 


(13) 


2.3.3.  VDC  Association  to  DFIG 

We  will  have  the  following  system: 


YVA] 

Vb 

Vpc  [  2  -1  -l] 
—  -  -1  2  —1 

r5ii 

Si 

[Vc\ 

3  Li  _!  2  J 

DiJ 

(14) 


In  our  case,  the  control  of  the  switches  of  the  inverter  is  performed  by  use  of  the  command 
modulation  or  PWM  pulse  width. 


3.  MAXIMUM  POWER  EXTRACTION  TECHNIQUE 

The  goal  of  the  (MPPT)  strategy  is  to  pick  up  the  maximum  power  from  the  wind;  it  involves  the 
following  of  the  power  curve  shown  in  Figure  4,  given  by  in  equation  (15): 

p„„  =  \c';r  (Xpt  )PA  V/,3  (15) 
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Figure  4.  MPPT  and  power  characteristics  in  function  of  mechanical  speed 


An  erroneous  speed  measurement  therefore  inevitably  leads  to  degradation  of  the  power  captured  by 
the  first  extraction  technique.  This  is  why  most  wind  turbines  are  controlled  without  control  of  the  speed. 
This  second  control  structure  based  on  the  assumption  that  the  wind  speed  varies  very  little  steady  state. 


APPLICATION  VECTOR  CONTROL  WITH  THE  ORIENTATION  FLOW  STATOR 

We  consider  the  DFIG  works  in  hypersynchronous  mode,  the  principle  is  to  direct  the  stator  flux 
along  the  axis  of  the  rotating  frame  [9] . 

So  we  have:  cpds  =  cps  and  we  have:  cpqs  =  0  Relations  between  current  stator  and  rotor  currents: 
The  technical  guidance  of  stator  flux  is  applying  on  the  couple  become: 


Ce  =  —2Pf(<Psd-Irq) 

If  neglecting  the  resistance  of  the  stator  winding  R_S,  voltage  expressions  become: 
ce  =  -3-pL-r(<Psd-irq) 


If  neglecting  the  resistance  of  the  stator  winding  R_S,  voltage  expressions  become: 


^sd 


0 


Ysq  —  Vs  —  (a)s.  <Psd 


[  7  _  (Ps  Lm  r 

L<zri  .  .  '  r 


(16) 

(17) 

(18) 

(19) 


4*s  —  ^sd  —  Ls-  Isd  Lm-  Ird 
0  —  Ls.  ISq  +  Lm.  Irq 

From  this  equation,  we  can  write  the  equations  linking  the  stator  currents  to  the  rotor  currents: 

(20) 

Relations  between  the  stator  and  rotor  currents  powers: 

[^s  ~  Ksd  ■  hd  Vsq  ■  hq 

Qs  ~  Vsq- hd  ~  Vsd-^sq 


/  = 

1sq  r  1rq 


(21) 


The  adaptation  of  these  equations  axis  system  chosen  and  the  simplifying  assumptions  made  in  this 
case  (Vsd  =  0  et  Vs  =  Vsq)  given: 


(Ps  =  Vs.Isq 

Ws  =  Is  -  hd 


(22) 
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By  replacing  the  stator  currents  by  their  values  we  get  the  expressions: 


Relations  between  rotor  and  rotor  currents  tensions: 

One  could  express  the  rotor  voltages  depending  on  rotor  currents,  we  can  write: 

Vrd  =  44  +  (Lr  ^ lrd  ~  9COs  (4  -  lrq 
vrq  =  Rrirq  +  (4  "4 )  ~at  ^rq  +  gc°s  ( Lr  ~  T f) lrd  +  R 

In  steady  state,  the  terms  involving  derivatives  disappear,  we  can  write: 

Vrd  —  ^r^rd  —  ^rq 

Vrq  =  RrIrq  +  goos  (4  - 17)  4  + 

Vector  control  in  stator  wind  power  system  based  on  a  PI  controller  [13,14]: 


(23) 


(24) 


(25) 


4.1.  Direct  Vector  Control 

This  control  mode  ensures  proper  decoupling  between  flux  and  torque  but  the  problem  this 
disadvantage  by  the  sensor.  We  neglect  the  terms  of  coupling  between  the  two  axes  of  control  because  of  the 
low  value  of  the  slip.  We  are  getting  as  shown  in  Figure  5. 
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Figure  5.  Block  diagram  of  the  direct  control 


4.2.  Vector  control  indirectly 

The  flow  is  controlled  in  open  loop.  It  is  not  measured  or  estimated.  The  quantities  (voltage  or 
current)  ensuring  the  flow  direction  and  are  evaluated  "decoupling  from  the  equations  of  the  machine 
transient. 

The  values  of  the  rotor  voltages  depending  on  power  and  found  we  calculate: 


V, 


dr 


v  =  - 


( j  Lm\ 

/  R  +(l  -dlL\ 

\ 

V  r  LS  )  n 

/  Rr+\Lr  Ls  J 

*'?  +  (; 

VsLm  r 

l  VsLm 

Ls 

V  LS 

/ 

l  l  -dn.\ 

Rr+\Lr  Ls)  r) 

*  P  gOJs' 

\Lr  os)*, 

VsLm  ° 

VsLm  *  ( 

Ls  J 

Ls 

\a)sLm 


V  Ls  J  0)sLm  J 


*Q  +  ga>s  (4  - 


0)sLm 


(26) 


The  Figure  6  shown  the  isolation  carried  at  the  outputs  of  regulators  in  rotor  currents  without  any 
return  to  the  system. 
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Figure  6.  Iindirect  command  without  power  loop  Figure  7.  Iindirect  control  with  power  loop 


The  Figure  7  shown  the  method  of  control  with  power.  In  this  method,  the  isolation  is  carried  out  at 
the  outputs  of  the  current  regulators  rotor  with  a  return  of  the  system. 


5.  RESULTS 

To  evaluate  and  test  the  indirect  control  technique  with  power  loop  of  active  and  reactive  power  by 
the  PI  controllers,  a  simulation  was  performed  under  the  MATLAB  /  Simulink.  In  this  case,  the  gains  of  the 
PI  controllers  are  based  to  method  of  design,  which  is  based  on  the  compensation  of  the  time  constant  of  the 
regulator  with  that  of  the  process  of  the  quantity  to  be  regulated,  and  were  refined  after  simulation: 

a.  For  the  power  loop:  k_p  =  75.75;  K  i  =  5354.55. 

b.  For  the  current  loop:  k_p  =  75.75;  K_i  =  5354.55. 


S  2  52  2  54  2  56  2  58  2  6  2  62  2  64  2  66  2  68  2.7 


Time(S) 


Time  (S) 


Figure  8.  Stator  current 


Figure  9.  Rotor  current 


Figure  10.  Speed  wind 


Figure  1 1 .  Rotor  Speed 
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Figure  12.  Stator  active  power 


Figure  13.  Stator  reactive  power 


Figure  14.  Comparison  between  Wr  and  Wr  with  MPPT 


6.  DISCUSSION 

We  obtained  good  results  in  dynamic  handling  and  response  to  regulations  imposed  and  reactive 
power.  The  fluctuations  in  the  power  due  to  the  PWM  inverter  and  the  dependence  of  these  powers  slip.  The 
figure  7  and  8  shown  the  results  of  the  stator  current  and  rotor  current  that  have  the  same  paces  with  the  trend 
of  the  wind  and  the  power.  Are  sinusoidal,  implying  a  clean  energy  without  harmonics  supplied  or  absorbed 
by  the  DFIG. 

The  results  obtained  in  Figure  12,  13  and  14,  we  can  conclude  that  this  type  of  control  is  more 
efficient  than  direct  control  in  terms  of  a  variable  speed  operation,  since  it  is  able  to  provide  a  decoupled 
control  of  active  power  and  responsive  regardless  of  the  drive  speed. 


7.  CONCLUSION 

This  paper  present  the  modelling  of  the  various  components  of  a  wind  system  for  distributed 
generation  of  electricity  and  the  study  of  different  double-fed  asynchronous  machine  control  systems  (DEIG) 
representing  generator  for  production  this  energy.  In  the  second  part,  we  have  begun  the  maximum  power 
extraction  technique  in  the  operation  of  the  wind;  this  method  proves  and  gives  good  results  for  the 
maximum  generated  power  to  the  grid.  Subsequently,  we  developed  vector  of  control  reactive  power  in  the 
stator  level,  the  proper  follow  instructions  for  the  two  powers  by  the  real  powers  debited  by  the  stator  of  the 
machine  showed  the  effectiveness  of  the  applied  control. 


NUMENCLATURE 

DFIG  Doubly  Fed  Induction  Generator 

p  Air  density 
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s 

turbine  area 

V 

Wind  Speed 

PC 

power  extraction  coefficient 

s(R) 

index  of  the  stator  (rotor) 

d,q 

indices  Park  repository 

V(I) 

voltage  (current) 

P(Q) 

active  power  (reactive) 

<D 

Magnetic  flux 

Cem  (Cr)electromagnetic  torque  (mechanical) 

R 

Resistance 

L  (M) 

Inductance  (mutual) 

a 

Coefficient  leaks,  a  =  1  -  M  2  /  LsLr 

Or  (Os) 

Position  of  the  rotor  (stator) 

cor  (cos) 

electric  speed  rotor  (stator) 

Q 

Mechanical  speed 

g 

Slip 

f 

Friction 

J 

Inertia 

P 

Number  of  pole  pairs 

Parameters  of  system: 

R  =35.25  m;  G=  90  m;  J=0.1Kg.m2. 
Vs=220/381;  Vr=  18/31  V;  f=50  Hz;  fr=14  Hz. 
R=0.3  ft;  L=0.014  H;  C=  2  e-3  F; 

Udc=620  V;  Rch=  85  ft;  E=220  V. 
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